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news and views 


ITH this issue, the TRANSACTIONS of the PGAP 

formally establishes itself as the organized tech- 

nical publication of the membership. The tech- 
nical papers that will be presented in this quarterly 
journal will have undergone the same careful and dis- 
criminating attention customarily accorded to papers 
published in the PRocEEDINGs. With the wholehearted 
support and cooperation of the Institute, we expect to 
see this new quarterly take its place as the official tech- 
nical organ for the workers in the Antenna and Propa- 
gation field. 

This section—News and Views—of the TRANSACTIONS 
will be a regular feature in each issue. Here, too, we hope 
to work out the problems of the Group by free discus- 
sion, welcoming free criticism as well as praise. We wish 
to use this section, in part, as a medium for ironing out 
controversial issues, a forum where the membership may 
feel free to express itself. 

A matter that we should like to bring to the attention 
of the Group covers the whole general subject of the 
form and organization of the professional groups—how 
many should there be, and how broad a field should each 
cover. The present organization reflects the thinking of 
arelatively small handful of energetic people who started 
with the basic premise that the Institute per se had 
grown much too big to satisfy the specialized needs of 
blocks of engineers. There is little doubt as to the wis- 
dom of this basic premise. However, are the professional 
groups, as currently constituted, properly and cohe- 
sively subdivided ? 

Listed below are the twenty-one professional groups 
recognized currently by the Institute. Truly, the IRE is 
going through a period of transition. 
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Airborne Electronics 

Antennas & Propagation 

Audio 

Broadcast & Television 

Receivers 

Broadcast 
Systems 

Circuit Theory 


Transmission 


_ Communications Systems 
~ Component Parts 
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_ Electronic Computers 


Electron Devices 


Engineering Management 
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Industrial Electronics 

Information Theory 

Instrumentation 

Medical Electronics 

Microwave Theory & 
Techniques 

Nuclear Science 

Quality Control 

Radio Telemetry & Re- 
mote Control 

Vehicular Communication 

Ultrasonics 


As an example of one opinion, we have received the 
following letter from Mr. Harold Wheeler of Wheeler 
Laboratories. His letter is reproduced here just as re- 
ceived: : 


“The time-honored association of Antennas and 
Propagation in IRE Committees has now outlived its 
usefulness. This has become more and more apparent 
when we note the striking division of interests within the 
Professional Group on Antennas and Propagation. 


“Therefore I propose that this Professional Group be 
split into two separate groups to be designated as fol- 
lows: 


Professional Group on Radio Antennas 
Professional Group on Radio Propagation 


“T think the need for such a separation will be apparent 
to anyone who is familiar with the group and thought- 
fully considers the question. If you should find a di- 
versity of opinion on this question of separation, then 
I should like to have an opportunity to argue the point. 


“Tt would be most beneficial if this separation could be 
made as soon as possible.” 


Here is a point of view that argues that the current 
subdivisions are not specialized enough. On the other 
hand, we are sure that there are many of the member- 
ship who feel that the present subdivision is already too 
fine. 


As one result of the formation of twenty-one separate 
professional groups, we now find a great increase in the 
number of technical meetings and conventions; more pa- 
pers are being presented and published, with wider geo- 
graphical distribution than ever before. The availability 
of organizing personnel in each professional group for 
sponsoring such meetings has provided this enlarged ac- 
tivity. But as with so many things, this may have its 
vices as well as its virtues. 


While it may well be that greater subdivision results 
in greater progress in each specialized area, it may 
equally well be that in the effort to nurture each further 
subdivision, the resulting wheat becomes increasingly 
indiscernible from the chaff. Then there is the problem 
of broad coordination of these various specialities. These 
are the questions vital to all of us in the IRE. 


We welcome your views for frank discussion. 
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ORGANIZATION 


PGAP Administrative Committee 


The operation of the Professional Group on Antennas 
and Propagation is under the guidance of the Adminis- 
trative Committee, as follows: 


Chairman: 
Dr. A. H. Waynick 
The Pennsylvania State College 
State College, Pennsylvania 
Vice-Chairman: 
MiP ey Carter 
RCA Laboratories Division 
Recky Points Wy NEY? 
Secretary-Treasurer: 
Mr. Delmer C. Ports 
Jansky and Bailey 
Washington 7, D. C. 
Administrative Committee Members: 
Dr. George Sinclair 
Sinclair Radio Laboratories, Ltd. 
Toronto, Canada 
Deni pony cl 
U.S. Navy Electronics Laboratory 
San Diego 52, Calif. 
Dr. H. G. Booker 
School of Electrical Engineering 
Cornell University 
Ithaca, N.Y. 
Mr. H. M. Wells 
Carnegie Institute of Washington 
Washington 15, D.C. 
Dr Leer van Atta 
Research and Development Laboratories 
Hughes Aircraft Company 
Culver City, Calif. 
Dre LafeChu 
Department of Electrical Engineering 
Massachusetts Institute of Technology 
Cambridge 39, Mass. 
Ex Officio Members: 
H. A. Finke 
Polytechnic Research & Development Co. 
Brooklyn 1, N. Y. 
J. S. Brown 
Andrew Corporation 
Sic hicaco 49, “LiL: 
Robert B. Jacques 
Electronics Division 
Thompson Products, Inc. 
Cleveland 14, Ohio 


PGAP Officers for 1953-54 


Mr. P. C. Carter has been elected as Chairman and 
Mr. D. C. Ports as Vice-Chairman for the year begin- 
ning June 1, 1953. The position of Secretary-Treasurer 
is open for the moment but will be filled very shortly. 


Membership 


All IRE members whose interest lies in the field of 
Antennas and Propagation are eligible for membership 
in the PGAP and will receive the TRANSACTIONS and all 
other group publications upon payment of the current 
assessment of $4.00 per year. The assessment is the 
Group’s principal means of covering the cost of the 
TRANSACTIONS. The TRANSACTIONS is sent to paid up 
members only. A glance at the present, initial issue 
should convince unpaid members that they would’ be 
missing a valuable contribution to their field of interest. 


Transactions 


The aim of the TRANSACTIONS is to serve as a profes- 
sional organ to meet and integrate the needs of the mem- 
bership. 

Contributed papers are given the same careful exami- 
nation as contributions to the PROCEEDINGS receive. 
The more restricted area of the TRANSACTIONS permit 
a prompt consideration of all papers. 

Dr. E. Weber, Chairman of the IRE Awards Com- 
mittee has indicated to us that papers which appear in 
Professional Group Transactions will be considered 
along with papers appearing in the PROCEEDINGS OF THE 
I.R.E. in connection with awards, such as the B. J. 
Thompson Award. The Executive Committee of the 
IRE also is considering setting up a policy wherein a 
$100.00 award would be given for the best paper pre- 
sented each year in each Professional Groups’ Transac- 
tions. 

It is hoped that the membership will take an active 
part in the success of our journal. Contributions are in- 
vited. , 

The diagram on the following page shows the admin- 
istrative organization of the TRANSACTIONS. 


GROUP CHAPTERS AND CHAPTER NEWS 


At the “grass roots” level, the professional group 
chapters should be the key to maintaining an active in- 
terest in the PGAP asa whole. In certain sections of the 
IRE where great interest in antennas and propagation 
is concerned, group chapters have been formed, are be- 
ing formed, or their possibility considered. PGAP mem- 
bers often, while realizing the value of such “grass roots” 
meetings and technical sessions, fail to lend their ener- 
gies in making local group meetings worthwhile and 
productive. 

Chapters have been established in Chicago and Los 
Angeles. Other chapters in the process of being formed 
are in the New York area, San Diego, and San Fran- 
cisco. PGAP members are urged to enter the activities 
of these chapters, and in areas where chapter activities — 
have not yet reached an organized stage, to investigate 
the ways and means for getting one going. J. Se Brown, 
Group Membership Chairman, stands ready to lend an 
accomplished hand in this matter. 
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Los Angeles. At the April 17 meeting of the chapter, 
attended by 62 persons, the following IRE Antenna and 
Propagation officers were elected: 


Chairman: 


Vice-Chairman: 


Secretary: 


Dr. M. J. Ehrlich 

Research Development Labs. 
Hughes Aircraft Co. 

Culver City, Calif. 

Dr. R. S. Wehner 

Research and Development Labs. 
Hughes Aircraft Co. 

Culver City, Calif. 

Mr. Robert Krausz 
Microwave Engineering Co. 
4000 Mount Lee Drive 

Los Angeles 28, Calif. 


- ' 
At this same meeting, Dr. John V. N. Granger of 
_ Stanford Research Institute gave a general talk on “Air- 
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Chicago. The following chapter officers were elected 
at the February meeting: 

Chairman: G. Kearse 

American Phenolic Corporation 
Vice-Chairman: L. R. Krahe 

Andrew Corporation 
Secretary: R. A. Jensen 

Sears Roebuck & Co. 

They will take office on July 1. 


The following chapter meetings and talks have been 
held during the past season: 


Oct. 17 —“Antenna Design for Microwave Relays,” 
by E. Dyk, Motorola, Inc. 

Nov. 21—“UHF Measurements and Smith Chart,” 
by Joseph Markin, Raytheon. 

Feb. 20—“The Corner Reflector Antenna,” by 
E. Harris, Mark Products Co. , 

Apr. 17 —“Antenna Problems in Air Navigation,” by 
S. Rachie, C. A. C. 


MEETINGS 


The group has participated directly or indirectly in 
several very successful meetings during the year and is 
now assisting in the arrangements for several more. dur- 
ing the forthcoming year. Dr. L. C. Van Atta is serving 
as Chairman of the Group meetings Committee. 

Future meetings of interest are: 


Conference on Radio Meteorology 


A conference on Radio Meteorology will be held at 
The University of Texas, Austin, Texas, from Novem- 
ber 9 through 12, 1953. Sponsoring and participating or- 
ganizations will include: 


American Meteorological Society—126th National 
Meeting. 

Radar Weather Conference—4th Annual Meeting. 

Professional Group on Antennas and Propagation of the 
Institute of Radio Engineers. 

National Commission II on Tropospheric Radio Propa- 
gation of the International Scientific Radio Union and 
the Joint Commission on Radio Meteorology. 


Papers are invited. It is tentatively planned to sched- 
ule review papers in order to permit more effective co- 
operation of the meteorologists and radio engineers on 
common problems in the field of radio meteorology. 

All communications and papers intended for the Con- 
ference should be sent to: 

Mr. John R. Gerhardt 
Conference on Radio Meteorology 
Box F, University Station 

Austin, Texas. 


Western Electronic Show and Convention, 
San Francisco, Calif. 


There will be PGAP sponsored sessions at this meet- 
ing to be held August 19-21, 1953. 


U.S. A. NATIONAL COMMITTEE 
URSI Sprine MEEtING, 1953 


(In collaboration with the I.R.E. Professional Group on 
Antennas and Propagation) 

An attendance of approximately 350 was obtained at 

the Spring 1953 meeting of the U.S.A. National Com- 

‘mittee of URSI, in collaboration with the I.R.E. Pro- 
fessional Group on Antennas and Propagation, held at 
the National Bureau of Standards, Washington, D. C. 
on April 27, 28, 29, and 30, 1953. 

The U.S.A. National Committee met on April 30, 
1953, and elected Officers for the years 1953-1955. The 
new officers are listed: 

Chairman: Dr. A. H. Waynick 

Pennsylvania State College 
Mr. H. W. Wells 

Carnegie Inst. of Washington 
Dr. W. E. Gordon 

Cornell University 


Vice-Chairman: 


Secretary-Treasurer: 
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Commission Chairmen: 
Commission 1 —Radio Measurements Standards: 
Dr. R. G. Fellers 
Naval Research Laboratory 
Commission 2 —Tropospheric Radio Propagation: 
Dr. H. G. Booker 
Cornell University 
Commission 3 —Ionospheric Radio Propagation: "he 
Dr. B. L. Berkner 
Associated Universities, Inc. 
Commission 4 —Terrestrial Radio Noise: 7a 
Mr. F. H. Dickson 
U. S. Department of the Nagy 
Commission 5 —Radio Astronomy: Qiring 
Dr. John P. Hagan 
Naval Research Laboratory 
Commission 6a—Antennas and Wave Guides: 
Dr J. Vi N. Granger 
Stanford Research Institute 
Commission 6b—Radio Waves and Circuits, Includ- 
ing General Theory: 
Dr. J. Pettit 
Stanford University 
Commission 7 —Electronics: 
Dr. J. R. Whinnery 
University of California 


A total of 83 papers were presented. Abstracts are 
available from the U. S. Secretary, Dr. W. E. Gordon, 
Cornell University, Ithaca, New York. No papers con- 
cerning the work of Commission VII, Electronics, were 
presented. Considerable interest was noted in the work 
of Commission V, Radio Astronomy, in which several 
noteworthy advances in this fiéld were reported. 

The next meeting will be held in the Fall of 1953 at 
Ottawa, Canada, at the invitation of, and in collabora- 
tion with, the Canadian National Committee of URSI. 
It is anticipated that in addition to the Professional 
Groups on Antennas and Propagation, other I.R.E. Pro- 
fessional Groups will also participate in this meeting. 

A listing of the commissions and papers which were 
presented under their sponsorship follows: 


Combined Session of Participating U.S.A. 
National Commissions 

Remarks concerning the Xth General Assembly, of 
URSI, Sydney, Australia, August, 1952—C. R. Bur- 
rows, Chairman, U.S.A. National Committee of URSE 
Cornell University, Ithaca, New York. 

Auroral Research at the Geophysical Institute, Uni- 
versity of Alaska—C. T. Elvey, University of Alaska, 
College, Alaska. 

The Probability Distribution of Atmospheric Noise— 
A. W. Sullivan, J. M. Barney and S. P. Hersperger, 
University of Florida, Gainesville, Florida. 

Interstellar 21 CM Radiation—H. I. Ewen, Harvard 
University, Cambridge, Massachusetts. 

Effects of the Moon on the Outer Atmosphere—A. G. 
McNish, Nat'l. Bureau of Standards, Washington, D. C. 
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contributions 


Measurement of Path Loss Between Miami 
and Key West at 3675 MC 
R. L. ROBBINS* 


Summary—Radio transmission measurements have been made 
at 3675 megacycles on the 130-mile path between Miami and Key 
West, Florida, which is largely over watery wastes of the Everglades 
and shallow sea waters of the Florida Keys. Path loss and fading 
characteristics for this terrain were not found to differ materially 
from the characteristics of hilly or mountainous paths in the north- 
eastern section of the country. 


HE PROPAGATION of microwave frequencies 
ike distances far beyond the optical horizon has 

been the subject of considerable study and ex- 
perimentation in recent years. K. Bullington! has sum- 
marized the principal findings published and has added 
data obtained by the Bell Telephone Laboratories over 
the rolling and mountainous terrain found in the North 
Atlantic and New England states. During September 
1952, the Laboratories conducted similar tests of 3675- 
mc propagation over a 130-mile path between Miami 
and Key West, Florida, where the transmission path 
was largely over the watery wastes of the Everglades 
and the shallow Gulf of Mexico waters north of the Flor- 
ida Keys. Somewhat less than fifteen per cent of the 
path was over dry land. 

Fig. 1 is a mass plot of all the available data on the 
median signal levels measured at various distances from 
the transmitter. The triangles represent previous find- 
ings and the dots represent data taken in Florida. 

The tests were conducted over five transmission paths, 
varying in length from 123 to 138 miles, between three 
transmitting sites in the vicinity of Miami and three re- 
ceiving sites at Key West. The sites chosen provided 
various foregrounds; i.e., water, relatively light vegeta- 
tion over smooth ground and the roof tops of residences. 


* Bell Telephone Labs., Inc., New York, N. Y. ' ; 
1K, Bullington, “Radio transmission beyond the horizon in the 


" 40- to 4000-mc band,” Proc. I.R.E., pp. 132-135; January, 1953. 


The first transmitting site, on the property of the Tropi- 
cal Radio Telegraph Co. about 13 miles north of down- 
town Miami, provided a clear transmission path over 
smooth ground having low vegetation for a distance of 
about one-half mile. The second transmitting site was 
located on the Rickenbacker Causeway adjacent to the 
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Fig. 1—Mass plot of SHF propagation data. 


island of Virginia Key. The smooth waters of Biscayne 
Bay extended southward from the very base of the tower 
for eight miles. The third transmitting site, at the west- 
ern border of Coral Gables, was on the roof of the “87” 
Exchange Building of the Southern Bell Telephone and 
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Telegraph Company where the transmission path was 
over the neighboring residences. The antenna was about 
forty-five feet above the ground. 

The first Key West receiving site was at the western 
edge of the Municipal Golf Course. The smooth man- 
grove-covered rocky terrain sloped gently into the wa- 
ter. The second receiving site was located on Bertha St. 
near the Atlantic Ocean. Here the signals arrived over 
the smooth shallow waters of a former salt evaporating 
pond and a group of residences about one mile away. 
The antenna tower was erected on level coral fill about 
75 feet from the water. The third receiving site was on 
the roof of the La Concha Hotel. The building is six 
stories high and the roof is well above the surrounding 
business buildings and residences. 

Fair weather, except for small scattered thunder- 
storms prevailed at both Miami and Key West during 
the period of the test. It is believed that similar weather 
existed between the two cities. The storms varied from 
light to sharp intensity. 

The received signals were observed and recorded for 
a period of approximately forty-eight hours during each 
test. As ascertained by observation of the receiving an- 
tenna alignment and of a magnetic compass, the signals 
were approaching over a great circle path. A +4° azi- 
muth adjustment of the 57-inch paraboloid receiving 
antenna (32 db gain) was required to reduce the re- 
ceived signal strength by 3 db. Such an antenna has a 
nominal beam width of 3.5° between the half-power 
points. A measure of directivity in the vertical plane 
could not be taken, but the signal was observed to drop 
rapidly as the antenna was tilted below the horizon and 
to drop relatively slowly as it was tilted above the hori- 
zon. Qualitative observations of directivity at the trans- 
mitter indicated that similar azimuth and elevation ef- 
fects were present there. 

The mieans by which the test signals were transmit- 
ted, recegved and evaluated will be briefly described 
later. It uffices now to know that they were pulsed sig- 
nals of 1°5 microseconds duration with a repetition rate 
of 400 cycles per second transmitted at a peak power of 
500 kw into a 32-db antenna. The received pulses were 
evaluated by a peak-reading vacuum tube voltmeter 
and recorded by a strip-chart recorder connected in se- 
ries with the voltmeter deflection coil. The recorder had 
a‘useful range of about 10 db. The time constant of the 
voltmeter-recorder combination was about 1.5 seconds. 

The received pulses were displayed on an oscilloscope 
while they were being recorded. Almost no lengthening 
of the pulses was observed; but they varied rapidly in 
amplitude and in shape, indicating transmission over a 
multiplicity of paths. Many of the fluctuations were 
probably too fast to be caught by eye. 

A typical strip-chart recording of the pulses is shown 
in Fig. 2. The lower portion of the chart was recorded 
at a speed of 15 minutes per division, or 3 inches per 
hour. The upper portion was recorded at a speed of 15 
seconds per division, or 3 inches per minute, to expand 
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the time scale. A calibration scale has been superim- 
posed at the bottom of the chart. It is seen that the pulse 
amplitude, as recorded, varied approximately +2 db 
about a short-term average. 
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Fig. 2—Typical strip chart recording of received signals. 


The curves of Fig. 3 are plots of the data taken by 
visually averaging, at 15-minute intervals, the record- 
ings made on each of the five transmission paths. Inter- 
ruptions in the curves are due to equipment failures. 
There were no signal fade-outs during the period of the 
tests. 

The variation of received signal power with changes 
of receiving antenna elevation is shown in Fig. 4. The 
solid line represents data taken by day and the dashed 
line is for data taken by night. Fig. 5 shows the effect of 
changing the transmitting antenna elevation. These 
curves, together with data taken at the “87” Exchange 
Building and the La Concha Hotel sites, indicate that 
there is no pronounced advantage to be obtained by ele- 
vating the antennas to as much as 100 feet. However, 
the relatively constant power received with changes of 
antenna elevation does encourage construction of larger 
directional antennas. 
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Fig. 4—Variation of received signal power with changes of 
receiving antenna elevation. 


The test gear comprised two vehicles and the appara- 

tus therein and two transportable steel towers with 57- 
inch paraboloidal antennas. A conventional radar 
transmitter of 500-kw peak power and its engine-gener- 

- ator power supply were contained in a 23-ton truck. 
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Fig. 5—Variation of received signal power with changes of 
transmitting antenna elevation. 


The transmitting antenna was fed through wave guide. 
Each of the antennas was mounted on gimbals so that 
remotely controlled motors could direct it in azimuth 
and elevation. In addition, each antenna was supported 
by a carriage so that it could be elevated to any position 
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on the tower. The 3675-mc receiving circuits were en- 
closed within a metal housing attached to the antenna. 
The IF, video, and power supply equipment was in- 
stalled in a station wagon. 

A block diagram of the receiving equipment is shown 
in Fig. 6. The microwave receiver unit contained a kly- 
stron oscillator, attenuator and wavemeter for control 
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nal generator containing a pulsed klystron oscillator, a 
wave guide attenuator and a directional coupler. The 
high-level output of the directional coupler was deter- 
mined by means of a thermistor-type power meter. The 
low-level output of the directional coupler was further 
reduced by two wave guide attenuators until it could be 
directly substituted for the signal at the antenna input 
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Fig. 6—Block diagram of receiving station equipment. 


and evaluation of the oscillator output, a crystal con- 
verter by which the incoming signal was demodulated 
to 60 mc, and a 60-mc preamplifier. The amplifier out- 
put was conveyed down the tower by coaxial cable to the 
remainder of the test equipment in the station wagon. 

The IF attenuator, IF amplifier, crystal detector and 
video amplifier provided the necessary gain and control 
thereof to permit operation and calibration of the vac- 
uum tube voltmeter and its associated recorder. The re- 
ceived pulses were observed on the oscilloscope as pre- 
viously mentioned. The equipment was calibrated by 
substituting a pulsed signal of like amplitude, duration 
and repetition rate, and of. known power, for the incom- 
ing signal. The calibrating pulse was secured from a sig- 


CONCLUSIONS 


1. On several beyond-horizon paths essentially over 
water in semi-tropical areas, the path losses observed 
were in excellent agreement with those reported in other 
areas by Bullington. 


2. No differences were observed from various fore- 
ground terrains. Overland, overwater and well elevated 
overland take-offs were investigated. 


3. Height-loss runs showed uniform received fields 
over the area that would be occupied by a large high- 
gain antenna. 


4, The diurnal and short-term signal variations were 
similar to those observed elsewhere. 
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Radiation From a Vertical Electric Dipole 
Over a Stratified Ground* 
JAMES R. WAIT} 


Summary—Expressions for the radiation fields at low frequencies 
of a vertical electric dipole situated on a horizontally stratified 
ground are derived. It is indicated that the well-known numerical 
results for the homogeneous ground can also be employed for ground 
wave propagation over a plane conductor composed of any number of 
parallel layers by suitably defining an ‘effective numerical dis- 
tance.” 


INTRODUCTION 


N THE THEORETICAL study of radiation from 
if an antenna it is usually assumed that the ground 

may be represented as a homogeneous medium with 
a specified conductivity and dielectric constant. The 
general problem has been investigated by Sommerfeld,! 
Van der Pol,? Niessen,? Bremmer,? Norton® and others. 
Extensive charts have been prepared by Norton’ which 
are very convenient for the determination of the field 
intensity of the antenna at some distant point on the 
surface of the earth. 

In this note it is indicated that the theory for the 
homogeneous earth can be usefully extended to the case 
of a horizontally stratified earth of any number of_ho- 
mogeneous layers. Because of the approximations in the 
analysis, the results are valid only for frequencies where 
the displacement currents in the air are small compared 
to the combined conduction and displacement currents 
in the ground. 


THE FORMAL SOLUTION 


The ground is considered to be suitably represented 
by a number of layers each with a thickness h, conduc- 
tivity a, dielectric constant e, and permeability pw in 
MKS units. A subscript 2 is added to these quantities to 
denote the mth layer from the surface of the ground. The 

- bottom or Nth layer is taken to be of infinite thickness 
and therefore is equivalent to a semi-infinite region. 

The region above the ground is considered to be a 
semi-infinite insulating space with dielectric constant €9 
and permeability wo. The emitter is represented by a 
vertical electric dipole of infinitesimal length ds and car- 
ries a current equal to the real part of I exp. (iwt). The 
dipole is situated in the insulating space at height h 
above the ground plane. 


* Work carried out under Project No. D48-95-II-14. (After writ- 
ing this paper the author came upon a paper by J. Grosskopf entitled 
“Das Strahlungsfeld eines vertikalen Dipolsenders uber geschichte- 
tem Boden,” Hockfreg. u Elek. 60, pp. 136-141; 1942. Grosskopf ob- 
tained a solution for the two layer ground which is a special case of 

_ the general result given in this paper.) 
+ Radio Physics Lab., Defense Research Board, Ottawa, Canada. 
1A. N. Sommerfeld, Ann. der Phys., vol. 81, p. 1135; 1926. 

; 2B. Van der Pol, Physica, vol. 2, p. 843; 1935. 

: 3K. F. Niessen, Ann. der Phys., vol. 18, p. 893; 1933. 

_ 4H, Bremmer, Phil. Mag., vol. 25, p. 817; 1938. 

5K. A. Norton, Proc. IRE, vol. 25, p. 1203; 1937. 


It is convenient to introduce a cylindrical polar co- 
ordinate system (p, ¢, z) such that the ground plane co- 
incides with the plane z=0 and the emitter is located on 
the z axis at z=h. Since there is symmetry about the z 
axis the fields can be represented everywhere by a Hertz 
vector with a component only in the z direction and is 
denoted by z. The fields are then expressed as: 


7 Or 5 
= y — ‘ 
*  Apdz : 
7? Or 
Es =O; Hy => = (1) 
tw Op 
Qe? 
E,=—-—-yat+ 1, lil = Oy 
02? 
where 
7? = iopw — euw?. 


The function 7 satisfies the wave equation: 
(V? — 7?)r = 0. (2) 


A subscript m is again added to the quantities to denote 
the layer that the quantities pertain to. The Hertz func- 
tion + can then be represented as a superposition of 
suitable solutions of (2) and is given as: 


Tn =f [A,,(A)ew? + B,(d)e~%# |Jo(Ap) dr, (3) 
0 
where 
Uy = (N+ Yn?) t?, 


Following Sommerfeld, the primary excitation of the 
dipole can be written as follows: 


Top = Cea /7y = all Nto le”) “OT (Ap) dr, (4) 
0 


for s<h, where r:= [z—h)?+p?]"”? and 
c = Ids/Ariwe. 
The coefficient Ao(A) is then given by: 
Ao(A) = ¢—)“0\ uot. 
Since only outgoing waves can be allowed in the bottom 
semi-infinite region (i.e. 2 =) it is seen that: 
By(d) = 0. 

The boundary conditions require that the E, and Hy 
field components are continuous at the interfaces be- 
tween the layers. These in turn require that the quan- 
tities 07/0, and ym are continuous at the interfaces. 
The 2N unknown coefficients A, and B, can therefore be 


easily expressed in terms of the known coefficient Ao by 
solving the system of 2N linear simultaneous equations 
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derived from the boundary conditions. Omitting the de- 
tails of the algebra the final result for the Hertz function 
qo in the air is given by: 


T) = T0p —+- cf uo RA) e— “0+ J o(Ap) dr, (5) 
0 


where 
Ko = Ve 


ays eae 
Z.+ Ki tanh mh 

; ie + Z, tanh uae 

_ 23+ Ke tanh Ugh 

j K+ Z3 tanh ahs 


R(A) = wodte%" Bo(d) = 


Li ae 


Jag, 


and so on, until 
Ky + Ky-1 tanh uy—shw-1 


Ky_1+ Kw tanh uy—ihy-1 


Zy—-1 = Kw-1 


where 


Kn = Unlon + twen)—. 


The formal solution of the problem is thus completely 
specified. The reduction of this result to a useful form is 
now undertaken. 


Tue HoMOGENEOUS GROUND 
When the ground is homogeneous such that h, tends 
to infinity, the factor R is simply given by: 
yuo — Your 


R(A) = r ; 
yuo + Yo" 


Sommerfeld! has shown that for the case of the emitter 
on the ground, (i.e. #=0) that #1 can be replaced by 11 
in the integral if | y2|>| vo?| so that 


co hem“ Tg(Ap) dd 
oe f (6) 
0 (M4 02)? + 02/11 


He also showed that this integral could be expressed in 
terms of the error function complement with a complex 
argument for the case | yor | >1 and z/p<1. His result 
can be conveniently written as follows: 


wo = 266009 
where F is an attenuation factor given by: 
F = 1 — i(rpi)\?e-”? erfc (tp2!!?) 


where 


Yo%p 


ay" 


Vi2\" 
po = pi (: =f ~*). 
Yop 
The vertical electric field #, in the air at the surface of 


the ground (i.e. z=0) with the same approximations, is 
then given by: 


Pie 


and 


duwl dse—VP 


LS ee 
in 0 (7) 
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and the horizontal field E, on the surface of the ground 
is given by: 


E, = (70/71) (8) 
where 
Fy = 1 — iv/rpi e-” erfc (ips). 
The function Fy has been tabulated by Norton? for a 


range of the values of both magnitude and phase of the 
complex factor 1. 
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Fig. 1—The amplitude of the wave correction factor Q for a stratified 
ground composed of two homogeneous layers. 


THE STRATIFIED GROUND 


When the ground is stratified such that there are NV 
layers present, the integral expression for mo is a great 
deal more complicated than the corresponding integral 
for the homogeneous ground discussed in the preceding 
section. The approximate solution, however, can be ob- 
tained in a similar manner by replacing wi, pe, Ms, etc. 
by 1, Yo ‘Ys, etc. This will be justified if | y2/y0?l, 
| -y22/ve?|, | v3?/Y0?| , etc. are much greater than unity. 
The expression for 79 for =0, is now given by: 

he~"2T (Ap) dd 


rT) = 2 f See 
oe Cw ioe male fy) avail wilt ie Ga ies 


where ye can be called an “effective propagation con- 
stant” and is given by: 


(9) 


Ver iuo/Z 
where 
Zo + N1 tanh yihi 


A=m 
mi + Ze tanh yh, 


1953 


Z3 + Ne tanh Yoho 
" n2 + Zs; tanh yohe 


Zr= 


and so on, until 


nv + nn-1 tanh yy_ihy_1 


Zy-1 = nN-1 
nv-1 + nw tanh yy_ihy_1 


where 
tn iue/Vn. 


By complete analogy to the case of the homogeneous 
ground, the electric fields E, and £, in the air at z=0 for 
| vop| >1 are given by: 


iuwl dse—ror 


eer rie (10) 
2rp 

and 

E, = (y0o/v)E« (11) 
where 

F, = 1 — t(rp.)'/"e-™ erfc (ip,1!/?) 
and 
Yo°p 
hehe Be Qy 2 


The factor p, can be called the “effective numerical dis- 
tance.” 

The “wave tilt” W is usually defined as the complex 
ratio of the horizontal field E, to the vertical field E,, 
and it is conveniently written: 


W=W.O 
where W>» is the wave tilt for the homogeneous ground 
given by: 

Wo = 0/¥1 
and Q is a correction factor given by: 


QO = "1/Ye. 


The “effective numerical distance” p, can then be writ- 
ten in terms of the numerical distance p; for the homo- 
geneous ground as follows: 


pe = iO”. 


It can be seen that when |y14:|>>1 the factor Q ap- 

proaches unity and the ground behaves as a homogene- 

ous half-space with electrical constants oj, €1, and p. 
The “Wave Tilt Correction Factor” for a two layer 


ground (i.e. h2= ©) takes the following form: 
vit v2 tanh yik1 

QO = | Q| = —— 
yo + v1 tanh yiky 


When displacement currents in the ground are negligible 
compared to the conduction currents, such that o>. 
and o2>ew then 


v1 = ioyw and Y2 ~ tow. 
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The amplitude factor |Q| and the phase gq are then 
plotted in Figs. 1 and 2 as a function of the “relative 
depth” V, given by V = (o,uw)/2h; for various values of the 
conductivity ratiog:/o,. It is interesting to note that when 
V<i and o)>0; that the “effective propagation con- 
stant” y. is somewhat greater than the propagation 
constant 71, of the upper layer. The “effective numerical 
distance” p, is also somewhat less than the correspond- 
ing value pf; for the homogeneous ground. At large 
“numerical distances” such that p, and p; are greater 
than about ten, a simple expression is obtained for the 
attenuation factor by considering only the first term of 
the asymptotic expansion of the error function comple- 
ment, so that: 


F,=—1/2p, and F,& — 1/2p, 
and then 


fost 'n/ O02. 


This implies that for this case the field intensity of the 
ground wave is considerably greater for the stratified 
ground than for the homogeneous ground. 
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Fig. 2—The phase of the correction factor Q. 


For deeper strata such that V>1 and o2>0; the field 
intensity is less for the stratified ground than for the 
homogeneous ground. For V greater than about 3.0 the 
factor Q becomes very close to unity, and the field in- 
tensity of the ground wave is determined only by the 
electrical constants of the upper layer. For a frequency 
of 125 kilocycles per second with a conductivity oi, of 
the upper layer equal to 10? mhos per meter, the depth 
h, corresponding to V=3.0, is about 30 meters. 

It is then concluded that the propagation of the 
ground wave of a radio-frequency transmitter is not 
usually influenced by irregularities in the earth’s crust, 
at depths greater than one hundred feet. A similar con- 
clusion was arrived at in a previous paper® where a two- 


dimensional incident wave was considered. 
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A Two-Dimensional Microwave Luneberg Lens 
G. D. M. PEELER¢, MEMBER, IRE, AND D. H. ARCHER{ 


Summary—A two-dimensional microwave model of the Luneberg 
lens has been designed employing the TE1o mode. It consists of two 
36-inch diameter, almost-parallel, conducting plates with the space 
between plates filled with polystyrene. Its thickness varies with the 
normalized radius, r, to give the desired index of refraction 
n=,/2—r. Due to symmetry about the center, this lens maintains 
constant gain and beam shape as a feed is scanned over its circum- 
ference, while the side lobe level remains at least 18 db below peak 
power. Experimental patterns show good agreement with computed 
patterns in the two principal planes. 


INTRODUCTION 


that can scan a radiation beam over wide an- 

gles in space. For low scanning rates, the entire 
antenna assembly can be moved to position the beam in 
space. As the scanning rate increases, moving the whole 
structure becomes impractical because of the large mass, 
and thus other methods must be used. One method 
widely used for obtaining higher scanning rates involves 
moving a low-mass feed over the focal surface of an ob- 
jective which is capable of maintaining good beam char- 
acteristics for various feed positions. The problem, then, 
is to obtain such an objective. 

The parabolic reflector, standard objective for non- 
scanning antennas, is geometrically perfect when fed at 
its focus and performs well within a few beamwidths of 
its axis; but its performance deteriorates rapidly for 
wider angles.! The spherical reflector, although not a 
perfect objective for any feed position, possesses more 
symmetry and can be used over wider angles when con- 
sidered with correctors such as the Schmidt correcting 
plate? or concentric-type correcting elements. The Lune- 
berg lens is an objective possessing complete symmetry. 
This lens is a sphere of unit radius whose index of re- 
fraction ” varies with the radius 7 according to the re- 
lation n=+/2—? as developed by Luneberg.* When fed 
by a source at a point on the surface of the sphere, this 
lens will focus a beam of rays at infinity (i.e., produce a 
plane wavefront). Since the focal surface of the lens is 
the surface of the sphere, the motion of a feed over the 
surface will produce a corresponding motion of the ra- 
diated beam in space without beam deterioration. 

This Luneberg sphere, used as a microwave lens to- 
gether with a suitable feed system, would in many re- 
spects be an ideal antenna for complex scans. However, 


VU ast RADAR SYSTEMS require antennas 


t+ Naval Research Laboratory, Washington, D. C. 
1K.S., Kelleher and H. P. Coleman, “Off-Axis Characteristics of 
me Es a Reflector,” NRL Report 4088 (Unclassified); Dec 
Ss 4 : 
2H. N. Chait, “A Microwave Schmidt System,” NRL 
3989 (Unclassified); May 14, 1952 i. eee 
3R. K. Luneberg, “Mathematical Theory of Optics,” 
Univ. Graduate School, Providence, R. I.; 1944, Ee 


such a three-dimensional lens cannot be easily con- 
structed at this time since there are no suitable mate- 
rials which will give the needed variation of index of re- 
fraction, although some of the expanded dielectrics and 
artificial dielectrics look encouraging. 

For many applications, it is not necessary that the 
antenna be capable of rapid scanning over a solid angle. 
An antenna capable of rapid scanning in a plane is quite 
often desirable, especially if the scan angle is large. This 
report is concerned with the development of an antenna 
to meet the latter need. 

Several models of the Luneberg lens which scan in a 


‘plane have been designed using parallel plate tech- 


niques. Rinehart! calculated a geodesic analogue which 
employs the TEM mode in a “derby-hat” shaped paral- 
lel-plate region. Experimental models using Rinehart’s 
design have been reported by Fine® at the Air Force 
Cambridge Research Center and by Kunz ef al.* Warren 
and Pinnell’ have designed the “tin hat,” a variation of 
the derby-hat which includes a bend in the mean sur- 
face near the edge of the lens. Jones® designed a model 
which used the TE) mode in air-filled, almost-parallel 
plates in which the index of refraction varied with the 
radius as m=no0\/2—r*, where mo <1. Refraction at the 
lens aperture (7%)#1) was corrected by adding to the 
edge of the lens a parallel-plate region terminated in a 
linear aperture. This addition maintains symmetry 
about only the diameter rather. than the center, result- 
ing in beam deterioration and loss of gain off-axis. Jasik® 
of Airborne Instruments Laboratory has designed a 
spherical Luneberg lens. 


Lens DESIGN 


The Luneberg lens reported here was designed by us- 
ing the technique of parallel plates propagating the TEio 
mode. It was developed in an attempt to eliminate some 
of the problems encountered in the above-mentioned 
designs. In order to explain its design, let us begin by 
considering a dielectric-filled rectangular waveguide. If 


*R.F. Rinehart, “A solution of the problem of rapid scanning for 
radar antennae,” Jour. Appl. Phys., vol. 19, p. 860; 1948. 

°E. C. Fine, “The Metal-Plate Analogue of a Luneberg Lens,” 
p. 30 of “Proceedings of the Third Symposium on Scanning Anten- 
nas,” Naval Research Laboratory; November 29-30, 1950. 

: K. = Kunz, ae Be and J. D. Johannesen, “Final Re- 
port on Contract -099-ac-141,” Case Institute of : 
September 30, 1949. Te ee 

‘PF. G. R. Warren and S. E. A. Pinnell, “The Tin Hat Scanni 
Rye A ane ce SNe iS of the Tin Hat Scanning Aotennae 

ictor Co. Ltd., Montreal, Canada, Technical R. 7 

July tband Sepe aceon chnical Reports 6 and 7; 
.S. D. Jones, “A Wide-angle Microwave Radiator,” 

Part III, vol. 97, p. 255; July, 1950. EY Se ae 

9H. Jasik, “Spherical Luneberg Lens,” p. 89 of “Proceedings of 

a Ye ea on Scanning Antennas,” NRL Report 4000; 

uly lo, . 
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this waveguide is propagating only the TE1) mode, the 
index of refraction is 


oN 4/ : aA. 
nN * € (), (1) 
where \= free-space wavelength; \, = waveguide wave- 
length; a=H-plane dimension of the waveguide; and 
e’=relative dielectric constant. Since m does not depend 
on the E-plane dimension of the waveguide, the wave- 
guide may be replaced by a pair of parallel plates. 

In this design the plate spacing is varied with position 
in order to achieve the proper variation of n. The as- 
sumption is made that this slight deformation from 
parallelism has an insignificant effect on the wave veloc- 
ity at any point, i.e., the index of refraction at a given 
point depends only upon the spacing at that point. This 
technique is not truly rigorous since the actual configu- 
ration of the plates should be considered as a boundary 
value problem, the solution probably being expressed as 
a sum of waveguide modes. The assumption being made 
is that the magnitudes of modes other than the TEj) are 
small enough to be neglected in this discussion. The con- 
cept is not new, for somewhat the same assumption is 
made in deriving the various modes in a waveguide with 
walls of finite conductivity. 

By combining (1) with the desired index of refraction 
for a Luneberg lens, n = »/2—?”, the resulting plate spac- 
ing is given by 


IN 
a= ————_-- 
2fe6 —24+77 


(2) 


This relation can be readily used to determine the place 
spacing for experimental models after the wavelength 
and the dielectric are selected. Examination shows that 
for r=0, (2) is meaningless unless e’ > 2. Other limits on 
the dielectric constant may be found by specifying op- 
eration above cutoff for the TE19 mode and below cutoff 
for higher order TE modes. For operation above cutoff 
in the TE. mode, let \<2av/e’ in (2) for r=1 (where 
cutoff will first occur) and obtain 1<~/e’/(e’—1) or 
e’>1, which imposes no additional restriction on e’. For 
operation below cutoff in the TEzo mode, let A>ave' in 
(2) for r=0 (the point of greatest difficulty) and obtain 
2>-/e'/(e—2), which is satisfied if ’ > 2.67. 

For our experimental model, designed for a wave- 
length of 3.2 cm, polystyrene was picked as the dielec- 
tric because it was readily available in large slabs, even 
though its tabulated dielectric constant of 2.48 was less 
than the minimum of 2.67 specified for operation below 
cutoff for the TEs. mode. It was assumed that there was 
little possibility of exciting the TE20 mode and that it 
could propagate only in a small region in the center of 
the lens even if it were excited. Equation (2) was used 
for determining the plate spacing as a function of the 
radius. This model, together with a sketch of the cross 
section which greatly exaggerates the curvature of the 

plates, is shown in Fig. 1. An aluminum plate was used 
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(b) 
Fig. 1—Photograph (a) and cross-sectional sketch of lens (6). 


on each side of the polystyrene, the thickness of the up- 
per and lower plates at the circumference being 0.38 and 
0.92 inches, respectively. The polystyrene slab was ob- 
tained by cementing together two pieces 11836 
inches with Monsanto polystyrene cement #106. The 
polystyrene seam forms a diameter of the lens. 

For ease of construction one surface of the polystyrene 
was kept plane while the other was machined to the de- 
sired contour. With this method of construction the 
mean surface of the lens is no longer a plane and could 
cause phase distortion, but the deviation from a plane 
is slight. The diameter of 36 inches was picked as the 
largest convenient for obtaining and machining both the 
aluminum plates and the polystyrene. This choice gives 
a diameter-to-wavelength ratio of 28.6 which is large 
enough to produce a reasonable beamwidth. Flanges for 
positioning a feed were arbitrarily made to cover 90 de- 
grees of the circumference and are located at the rear of 
the lens. The half-ring at the front of the lens was added 
to the upper plate so that the dielectric was symmetri- 
cally located, making the total thickness 2.36 inches. A 
small amount of preliminary data indicated improper 
plate spacing due to a slight warp in the polystyrene 
slab which occurred after machining. To provide better 
control of plate spacing bolts were added a few at a time 
until the data showed that the warping was corrected. 
These bolts are perpendicular to the electric field and 
cause negligible reflections. 


LENS ILLUMINATION 


The feed for this lens is shown in Fig. 2. The a-dimen- 
sion of the waveguide aperture (0.52 inches) was chosen 
to fit the edge of the lens, and the 0-dimension to pro- 
vide a desirable lens illumination. The waveguide 
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aperture was filled with polystyrene to permit operation 
above cutoff. Both the waveguide and the polystyrene 
were tapered for impedance considerations, the length 
of the tapers being chosen to keep the variation of wave- 
guide wavelength as small as possible in this transition 
region. The two metal fins were added in the aperture to 
suppress the TEo: mode which would excite the TEM 
mode in the lens. The free-space E-plane radiation pat- 
tern for this feed is shown in Fig. 3. Cross polarization 
(presence of the TE: mode) is at least 25 db below peak. 


ES 
VN 


.90" (b) 


Fig. 2—Feed horn for lens. 


A relationship must be found between the distribution 
of the field from this feed and the field at the lens aper- 
ture in order to determine the illumination function at 


the aperture which will be used later in radiation pat- 


tern calculations. This can be done by first considering 
a typical ray in the lens as shown in Fig. 4. The point P, 


is defined as the point on the ray nearest the origin, O. 


Due to the symmetry of the lens, the path of the ray is 
symmetrical about P; so that the “triangles” OP)P, and 
OP2P; are congruent, and B=+y. Since the ray at P» is 
parallel to OP;, y=W=8. 
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Energy flow through the lens must now be consid- 
ered.!° From principles of geometrical optics it can be 
shown that energy flow is along rays, and if a tube of 
rays is considered, that the energy flow is constant 
through various cross sections normal to the tube. Then 
the energy flow through the wedge-shaped tube at the 
feed between 8 and B+d8 and between x and «+dx is 
Si(B, x)d8 dx, where S,(B, x) is the energy distribution 
per unit angle in the primary pattern of the feed, and x 
is the co-ordinate perpendicular to the plane of the lens. 


fsa} a 
a 

\ / 
z* * 
x 
> 10} 
(o) 
2 Lawl 

— E-PLANE FEED PATTERN MEASURED 

an IN FREE SPACE 
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LE t =( FEED PATTERN) .( SEC yw) ie x 
= —w— ASSUMED E-PLANE FEED PATTERN \ 
Lu 20 os as seas ae 
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Fig. 3—Feed pattern and aperture illumination. 


At the lens aperture in the region of P; the energy flow 
through the tube is Sa dy dx, where S, is the energy flow 
per unit area at the aperture. Since S, is the magnitude 
of the Poynting vector at the aperture, it will be as- 
sumed that it is proportional to the square of the mag- 
nitude of the electric field E,-at the aperture. In Fig. 4 
it is seen that y=ro sin y. Then, equating these two en- 
ergy flows (remembering that B=w), the relation be- 
tween the aperture field and the feed pattern is 


| Ea | = CrV/Si sec y, (3) 


where C; combines the multiplying constants. This ap- 
erture field is also plotted in Fig. 3. 


Fig. 4—Typical ray in lens. 


_ 3° The background for these considerations is di 

ia Siva Si Sia ( Mictpyaye Antenna Theory and Degas > ene 
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The problem of approximating S; with a function that 
lends itself to mathematical calculation will now be con- 
sidered. It will be assumed that S; is a separable func- 
tion of x and wy of the form 

Si(y, ¥) = S2(a)-S.(y). (4) 
Since the lens employs the TEi) mode, S2(x) is known 
to be 


S2(x) = Cz cos? 


TX 5 
ey (5) 


where ¢ is the spacing between plates at the lens aper- 
ture and C, is a constant. An expression for S3(wW) will be 
chosen to fit the feed pattern as closely as possible. As 
shown in Fig. 3, if 


Ss(v) = cos* y (6) 
this function fits the feed pattern quite closely for most 
values of y. It was not felt advisable to obtain a closer 
fit at the wider angles, since there is some doubt as to 
the exact value of the lens illumination in the region af- 
fected by this low-intensity field. The reason for ques- 
tioning the illumination in this region can be under- 
stood by considering the second curve of Fig. 3, which 
represents the feed pattern multiplied by the space at- 
tenuation factor, sec y, obtained from geometrical op- 
tics. Apparently geometrical optics is not valid in this 
region. It is felt that the field does not accumulate at 
the lens edges and then drop sharply to zero, but that 


fa 
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some of the energy predicted to appear at the edge of the 
lens actually appears outside the lens. Therefore, it is 
felt that the cos? y approximation to the pattern is suf- 
ficiently accurate. By substituting (4) through (6) into 
(3) and dropping constant multiplying factors, 


(7) 


| one = cos = cos y. 


RADIATION PATTERN CALCULATIONS 


Several of the Luneberg lens reports have considered 
the problem of calculating the radiation pattern. War- 
ren and Pinnell’ obtained patterns by numerical inte- 
gration. Although Walsh" avoided numerical integra- 
tion in a similar problem on cylindrical arrays by ap- 
proximating the array of discrete elements with a‘con- 
tinuous current sheet, it was felt that the Luneberg pat- 
terns could be obtained in an easier, more direct man- 
ner. 

Since rays emerging from the lens are not perpendicu- 
lar to the aperture, and consequently the field vectors 
are not parallel to the aperture, it is necessary to con- 
sider an approach using these fields for formulating ex- 
pressions for the radiated field rather than using stand- 
ard equations containing only scalar functions. The de- 
tails of these considerations and the calculations (see the 
Appendix) result in the following expressions for the 
normalized patterns: 


11 J. E. Walsh, “Radiation patterns of arrays on a reflecting cylin- 
der,” Proc. I.R. E. , vol. 39, p. 1074; 1951. 
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Fig. 5—E-plane calculated and typical experimental radiation patterns. 
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For the E-Plane: 


Ba(0) = cos? — [In(B) + (B)] = 5 


6 6 6 
-| sin B cos 5) — Boos os (B cos ko) 
Z 2 2 


where B= 2kro sin (0/2), R= 2a/d and 79 =lens radius. 
For the H-Plane: 


1 + cos ‘) (= sin *) 
——— } cos 
( 2 2 
kt sin 6\2 
1 = 
Tv 


a 
D 


E,(6) = 


[fan HO) 


where D = 2kry sin? (6/2) and t=aperture height. 

These radiation patterns, given above by (8) and (9), 
are plotted in Figs. 5 and 6, where they can be com- 
pared with experimental patterns. 
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Fig. 6—H-plane calculated and typical experimental 
radiation patterns. 


BEAMWIDTH FORMULAS 


The preceding solutions of the radiation field will now 
be used to develop formulas for the half-power oon 
widths of the Luneberg lens. 

The E-plane radiation pattern given in (8) is rather 

‘unwieldy for this purpose, but considerable simplifica- 
tion occurs if the range of values which 7o/A may assume 
is restricted. Calculation has shown that for r9/A=7 the 
effect of the second term in this equation is completely 
negligible and cos? (0/2) may be replaced by unity. 
Therefore, with 7o/\ sufficiently large, the expression for 
the radiation pattern becomes 


E(8) = Jo(B) + Jx(B) = f(B). 


A plot of E(@) versus B gives a universal radiation pat- 
tern valid for a lens of any size. From this pattern the 
beamwidth can be expressed in terms of B or 7o/X. How- 
ever, this same information can be obtained from the 
plot of (8) for the experimental lens by noting that the 
half-power points occur at angles 6;= +1.03°. Thus, for 
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any lens, 6:=2 sin! (B,/2kro), and the full beamwidth 


is given by 
r 
To 


in degrees. This equation is plotted in Fig. 7. 

Calculation has shown that for 7<7o/AS70 the first 
term in (9) may be approximated by a constant. Thus 
the H-plane radiation pattern may be written as E(@) 
=g(D), and a plot of E(6) versus D shows that the half- 
power points occur at angles #;= +19°. Thus, for any 
lens, 6,=2 sin-!4/D,/2kro, and the full beamwidth is 
given by 


(10) 


Gy 2 U3/\ra (11) 


in degrees. This equation is also plotted in Fig. 7. 


& 60% 


S 
fo} 
° 


E-PLANE BEAMWIDTH 
cy % cy 
H-PLANE BEAMWIDTH 
wo 
° 


‘0 0 20 30 40 50 60 70 
LENS RADIUS (IN WAVELENGTHS) 


Fig. 7—Beamwidths versus 79/A for E- and H-planes. 
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EXPERIMENTAL EVALUATION OF THE LENS 


Phase Measurements 

The phase front in the plane of the lens was measured 
using standard techniques.” A typical set of data is 
shown in Fig. 8. For a range of pickup positions of 25 
inches the maximum deviation from a linear phase front 
is 0.025 wavelengths. 
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Fig. 8—Typical phase front of lens. 


E-Plane Patierns 


Early E-plane patterns taken with the lens assembled 
without bolts had an unsymmetrical beam shape due to 
the slightly warped polystyrene which caused incorrect 


For example, see Silver, S., op. cit., sec. 15.11. 
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Fig. 9—E-plane radiation patterns for various feed positions. 


sis 
Sm (hune 
mse: 


rey 
al 


PO RR aL 


Ce ey 
Fr Grif Bn 


2 
Wi 


ERRGE am) 
SE (oad NB 


40° 
Sr eS OL 


Baca cate 
ies sta 
os es 


320 


330 


10 


20 


30 


Peeler and Archer: A Two-Dimensional Microwave Luneberg Lens 


Lf 


plate spacing. As mentioned previously, the addition of 
the bolts corrected this difficulty. Patterns taken with 
the polystyrene seam in various positions show no ap- 
preciable effects due to its presence. 

A typical E-plane pattern for the feed at zero degrees 
shown in Fig. 5 was a half-power beamwidth of 2.2 de- 
grees and side lobes 19.4 db below peak power. The 
beamwidth factor, defined as (beamwidth Xlens diam- 
eter) /wavelength, is 63 degrees for this lens, which indi- 
cates high efficiency in this plane. This pattern agrees 
quite well with the calculated pattern, plotted in Fig. 5. 

E-plane patterns for feed positions in increments of 
2.5 degrees around the circumference of the lens are 
shown in Fig. 9. Over the 90-degree range permitted by 
the flanges, the gain and beam shape remain reasonably 
constant, and most of the side lobes are at least 18 db 
below peak power. The angle of beam shift is equal to 
the angle of feed shift. 


H-Plane Patterns 


Fig. 6 contains plots of a typical H-plane pattern for 
the feed at zero degrees and the calculated pattern given 
by (9). These patterns have flat tops and no side lobes 
or minima as such, but the skirts extend to rather large 
angles. The half-power beam widths are 42.8 and 38 de- 
grees for the experimental and calculated patterns re- 
spectively. 

Because of the end-fire effect present in this plane, 
these beamwidths are much smaller than would be ex- 
pected from the H-plane aperture of 0.52 inch. This end- 
fire effect might be seen in the calculated pattern given 
in (9) if the first bracket is considered as the “element 
pattern” and the second bracket as the “array factor.” 
The H-plane beamwidth formula given by (11) shows 
that, due to this end-fire effect, the beamwidth varies 
inversely as the square root of the lens radius. Therefore, 
if the lens radius were increased to obtain a smaller E- 
plane beamwidth, the H-plane beamwidth also becomes 
smaller; however, for any practical radius the lens pro- 
duces a fan-shaped beam. If a more directive beam is 
desired, it is necessary to increase either the H-plane 
aperture or the aperture efficiency, or both, i.e., by ob- 
taining a more directive “element pattern.” 


Addition of Flares 

In order to increase the lens aperture, the addition of 
flares was considered, keeping in mind as requirements 
for the flares: symmetry about the center of the lens, 
freedom of feed movement, and noninterference with 
E-plane patterns. Since this lens uses the TE;) mode, the 
addition of any flares makes the index of refraction less 
than unity within the flares. Since Luneberg’s resulting 
variation of index of refraction inside the lens is based 
on the lens being imbedded in a medium with an index 
of unity, the effects of the addition of flares should be 
compensated for in a redesign of the lens. Although 
analytical attempts to include flares in a lens design 
were unsuccessful, the experimental addition of flares 
proved somewhat fruitful. 


18 L.R.E. TRANSACTIONS—ANTENNAS AND PROPAGATI ON 


Two sets of flares were investigated. A photograph of 
one set mounted on the lens is shown in Fig. 10. These 
“conical flares” are composed of plane sections cut to a 
radius of 20 inches (two inches larger than the lens ra- 
dius) to which are soldered 7-inch sections of a 90 degree 
cone, making the total H-plane aperture approximately 
12 inches. These flares increase both the aperture size 


Fig. 10—Lens with conical flares. 


and its efficiency, the greater efficiency resulting from 
the box horn! formed by the lens and the plane sections. 
This conical section could not be extended to the actual 
aperture of the lens since the spacing between flares 
would then be small enough to make that region below 
cutoff. Mounting slots permit adjustment of the flares 
relative to the lens. An H-plane pattern for the best po- 
sition of the flares is shown in Fig. 11, where it can be 
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Fig. 11—H-plane radiation patterns for flares in 
best position. 


13 See Silver, S., op. cit., sec. 10.15. 
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compared with other patterns. With these flares the 
half-power beamwidth of this pattern is decreased to 22 
degrees and side lobes are 18 db below peak power. 

A second set of “plane flares,” made identical to the 
plane sections of the conical flares and which provide an 
H-plane aperture of 2.36 inches, was also tried. An H- 
plane pattern with these flares adjusted to the best posi- 
tion is also shown in Fig. 11. The half-power beamwidth 
of this pattern is 24 degrees, and side lobes are 16 db be- 
low peak power. 

A comparison of patterns in Fig. 11 shows little dif- 
ference in the patterns taken with the two sets of flares. 
The larger aperture of the conical flares obtains little 
reduction in beamwidth or side lobe level and increases 
the gain by less than 1 db. It is felt that the smaller size 
of the plane flares will in most applications make them 
more useful. In the next section, concerned with stack- 
ing these lenses to obtain further reduction in H-plane 
beamwidth, these smaller flares must be used. 

As was pointed out earlier, the addition of flares alters 
the optics of the lens; however, E-plane patterns taken 
with the conical flares for various feed positions as 
shown in Fig. 12 indicate little variation from those in 
Fig. 9 for the lens without flares. 

Since these flares improve the H-plane patterns with- 
out appreciable effect on E-plane patterns, they can also 
be utilized as flanges for positioning the feed. Thus, by 
extending them over the full circumference of the lens, 
full 360-degree scanning can be accomplished. 


Stacked Lenses 


Since the use of flares to sharpen the beam of this lens 
proved limited, a further reduction in the H-plane beam- 
width was obtained by stacking several of these lenses 
one above the other and feeding them in phase. The 
stacking of two lenses will now be considered. Letd be the 
separation between centers of the apertures and 6 be the 
H-plane pattern angle as indicated in the Appendix. 
Then at the angle @ in the H-plane there is a phase dif- 
ference of (27d/h) sin @ between the fields from the two 
lenses. If E(6) is the field of a single lens, this combina- 
tion will produce the field 


E,(6) = E(0) + E()e(i274 sin 6) 


ad sin 0 


2E(O)e(ird sin 0) /X cos 


If E(@) is approximated by e—, where cz is a constant 
which is equal to 0.00347 for a single lens beamwidth of 
20 degrees, this expression becomes, when considering 
only the normalized magnitude, 
E,(0) = e~0.0034762 cos md sin 6 : 


¢ 


(12) 


. . ‘ y 
The resulting beamwidth and side-lobe level are thus a 


function of the separation, d. For any given value of d, 


the half-power beamwidth is obtained by equating 
E;(6) to 0.707 and solving the resulting transcendental 
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Fig. 12—E-plane radiation patterns for various feed positions with conical flares. 
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equation for 61, the half-beamwidth. The m** side lobe 
occurs approximately at 9=sin—! (md/d), and its level 
can be found by substituting this angle in (12). For this 
double-lens combination, the half-power beamwidth and 
the first side-lobe level for several values of d are given 
in Table 1. 


DA BLES 


H-Plane Half-Power Beamwidth and First Side-Lobe 
. Level for Double-Lens Combination 


d Half-Power Beam- First Side-Lobe 
width (degrees) Level (db) 
2r 11.8 26.9 
BEN Op 16.7 
3X 8.7 11.4 


Contour Patterns 


Since a phase front from this lens is not a plane, but 
somewhat saddle-shaped, it is not immediately apparent 
that the largest side lobes are found in the principal 
planes. For this reason, other sections through the ra- 
diation beam were explored. A contour plot of the beam 
structure from the lens without flares is shown in Fig. 
13. No side lobes appear which are greater than those 
found in the principal planes. Measurements taken for 
the lens with flares gave a similar result. 


Pattern Bandwidth 


All the data presented thus far were taken at a wave- 
length of 3.2 cm. However, Table 2 gives E-plane pat- 
tern characteristics for a range of frequencies from 8,900 
mc to 9,800 mc for the lens without flares. Within this 
9.6 per cent band the patterns are essentially constant. 
At lower frequencies, the beam shape begins to deterio- 
rate with an accompanying increase in the side-lobe 
level. Equipment was not available for taking measure- 
ments at higher frequencies. Similar results were ob- 
tained for the lens with flares. 


TABLE 2 
E-Plane Pattern Characteristics for Various Frequencies 
Frequency Beamwidth Side-Lobe eas 

(mc) (degrees) Level (db) (degtees) 
8900 3,1 16.9 88.6 
9000 2.8 Bale 80.0 
9100 26 ise) 71.4 
9200 Sher: 19.9 65.7 
9300 225 i faks 64.3 
9375 BS25 20.4 64.3 
9500 22 18.8 62.9 
9600 23 19.5 65.7 
9700 2.4 20.8 68.6 
9800 2.6 19 o7 74.3 


Impedance Match 


For the lens without flares, the vswr in the feed line is 
below 1.4 over most of the range of \ from 3.1 to 3.5 cm. 
No attempts were made to obtain an impedance match 
other than putting tapers in the feed as shown in Fig. 2. 
However, since the characteristic impedance of the lens 
at its edge is equal to that of free space at the design 
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Fig. 13—Radiation pattern contour plot (in 3 db steps) for lens with no flares. 


wavelength, the lens possesses an inherent match in this 
region. Mismatch between feed and lens was not inves- 
tigated except to decide that no choke is needed at this 
junction since currents for the TE19 mode are parallel to 
the gap between the plates of the lens and the feed. 


Gain 


The difference in peak power in the patterns of Fig. 
11 gives a comparison between the gain of the lens with 
and without flares. The measured gain at \=3.2 cm was 
22.2 db for the lens only and 25.9 db for the lens with the 
plane flares. 

There are several bases for evaluating the gain from 
pencil-beam antennas and, although this antenna does 
not produce a pencil beam, such methods will be used. 
The first method is to compare the gain with that ex- 
pected from a uniformly illuminated aperture which is 
47A/)*, where A is the aperture area. For this lens with 
the plane flares, the projected aperture is 402.36 
inches, so that at \=3.2 cm, 47A/\?=750=28.8 db. A 
comparison of the measured gain of 25.9 db =389 with 
this figure shows that the aperture efficiency is 0.52, 
which is good for a wide-angle scanner. 

The second method is to compare the gain with 
27,000/6z 84, where 9% and 6y are the E- and H-plane 
half-power beamwidths in degrees, respectively. If the 
beamwidths for this lens with the plane flares are sub- 


stituted in this expression, it becomes 27,000/(2.2 X24) 
=512=27.1 db. There is a difference of 1.2 db between 
this figure and the measured gain. 

If account is taken of a small loss in the polystyrene of 
the lens and also the power contained in side lobes which 
in the E-plane do not fall off rapidly with angle, this 
lens appears quite efficient. 


CONCLUSIONS 


A two-dimensional microwave analogue of the Lune- 
berg lens which utilizes the TE:9 mode has been de- 
signed which produces good radiation patterns and 
maintains good pattern characteristics over its limit of 
scan and over a band of frequencies. The calculated and 
experimental patterns agree quite well. 

This lens design lends itself to easy construction by 
requiring machining of only one surface of the dielectric, 
resulting in a lens with a small volume which is quite 
rugged. The weight of this lens can be reduced to that of 
the dielectric alone by coating the dielectric with a thin 
metallic layer and eliminating the heavy aluminum 
plates and the bolts used in the model described hereto- 
fore. The shape of this lens permits the addition of 
flares and the stacking of several lenses to produce a 
more directional beam than a single lens possesses. The 
flares may be used to position the feed, thereby permit- 
ting scanning over the full 360 degrees. 
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Fig. 14—Co-ordinate systems for radiation pattern calculations. 
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index of refraction of unity in the aperture, E, and H, 
Apeanpix are related as in a plane wave. Thus the aperture field 
may be expressed as 
RADIATION PATTERN CALCULATIONS 
The aperture-field method will be used to calculate EE, = — ty| Eq| ef#ro(i—eos ¥), (15) 


the far-field radiation pattern of the Luneberg lens. The 


solution appears in the literature,'4 and the components Pepe! fa Perr: 
of the radiation field are H, = t.| Ea| i eg . (16) 
fois = arenes From Fig. 14 it follows that 
Eg= — eine f | xt te = axbets| cr 
4a z ag p = t2% — tyro sin Y + 2.70 cosy, (17) 
“eile Mads, (13) R, = 7, sin 6 cos ¢ + ty sin @ sin @ + 7, cos @, (18) 
j <t A rte tp = 7, cos 0 cos ¢ + 7, cos O sin @ — i, sin 8, (19) 
Mo : eS in ~ ,, 
ae cite f | axe tot © 5xBa: | to = — i, sin o + 7, cos ¢, (20) 


nm = — t,sin y + 7, cos y. (21) 


a — If (15) through (21) are substituted into (13) and (14), 
where E, and H, are the values of the field in the aper- they take the form 
ture and 7 is the outward-drawn unit vector normal to 


-eike-Rids, (14) 


the aperture. Referring to Fig. 14, these expressionsmay 7, _ _ Jorov Ho€o en ikR f | Eq | ef#Rat—eos v) 
be simplified by the following considerations. Since the° 4rR j 
lens was designed to produce constant phase in the plane -{cos w sin #(cos 6 + 1) — siny sin o} eit Bidydx (22) 


of the lens along the line tangent to the aperture center, 


joror/ Moeo : : 
E, and H, have a phasing of kro(i—cos W) where Ey = — ae cine f | E.| eikro(1—cos ¥) 


4rR j 


“4S, Silver, op. cit., p. 161, equation (103). -cos wy cos (1 + cos dei? Bidydx. (23) 


, 
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Substituting (7) for [ee into these two equations and 
dropping constant multiplying factors, they become 


t/2 Tx w/2 ; 
hy= f cos —— ax f cos W {cos p sin ¢(1 + cos 4) 
us t a 


a /2 


— sin w sin 0} em ihre cos veike Ridy, (24) 


t/2 rx 1/2 
E, = ii Cospeeee ax f cos? y cos ¢(1 + cos @) 
at t os 


w/2 


-e7tkro cos Veikp ‘Rid, 


(25) 


Since patterns in the principal planes are of greatest in- 
terest, these special cases will now be considered. 
E-Plane Pattern 


For the E-plane pattern, the far-field point P is con- 
fined to the yz-plane where ¢= +90°, but due to sym- 
metry about this plane, only one value of ¢ need be con- 
sidered. Inspection of (17) and (18) for ¢ = — 90° shows 


> ~ 
oR, = ro sin W sin 6+ 79 cos W cos 8 = ro cos (y — 6). 


Then (24) and (25) may be written 


E4(@) 


t/2 rx w/2 
i cos =n av f cos [cos (W—8)-+cos y| 


t/2 w/2 
- et krolcos(y—8)—cos WI dy, (26) 


E, = 0. 


Since the integral in x in (26) is independent of 6, it isa 
constant as far as the radiation pattern is concerned, 
and thus may be neglected. The second integral may be 
rearranged to give 

r/2 


cos w cos (W — 6/2) 


-e2ikro sin 0/2 sin (Y—-9/2) dy, 


E,(6) = cos a/2 f 


—1/2 


(27) 


Making the substitutions B = 2kr) sin 6/2 anda=y—6/ 
2 gives 


a /2—0/2 fs) 
Eo(8) = cos 0/2 f cos (« -+ *) cos aeiB sin adgy, 
—1/2—8/2 2 
Q§ pe 7i2—8/2 
E4(0) = cos? a cos? aeiB sin «dg 
<r /2—6/2 
sin 0 a /2—0/2 
=e f sin a cos ae?B sin ada, 
2 —1r/2—6/2 


6 sin 6 
E,(0) = (cos? -) Ii — ( Yr (28) 
2 2 
The first integral, 1, is quite difficult to evaluate, but 
by dropping the 0/2 term from the limits of integration, 
it can be expressed as follows 


r/2 
t= f cos® we sin dq 


a [2 


w/2 
= f cos? a cos (B sin a)da 


a /2 
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w/2 

+ if cos? a sin (B sin a)da. 
—7/2 

The integrand in the second integral above is an odd 

function, so the value of the integral is zero. The inte- 

grand in the first integral is an even function, so it may 

be written as 


x/2 
I; 2 f cos? a cos (B sin a)da, 
0 


w/2 
Ih= f cos (B sin a)da 
0 
7/2 
+ ip cos 2a cos (B sin a)da, 
0 


1 =~ [IB) + J2(B)]. (29) 
For values of 6 up to 6°, which is sufficient to include the 
first two side lobes, the error introduced by neglecting 
9/2 in the limits of J, can be shown to be less than 0.004 
per cent of (I)max. The second integral, J2, in (28) can 
be integrated directly to give 


277. 6 
I, = —| sin {| B cos— 
5% 2, 
0 0 
— Bcos— cos{ Bcos—)|]. (30) 
2 2 


Substituting (29) and (30) into (28), the normalized 
radiation pattern for the E-plane given in the text as 
(8) is then 


£ 
4 


2 sin 0 
a B2 


8 4 8 
sin (2 cos -) — Bcos— cos (2 cos -)| . (8) 
2 2 2 


H-Plane Pattern 


For the H-plane pattern, the far-field point P is con- 
fined to the xz-plane where 6=0 or 180°, but due to 
symmetry about this plane only one value of @ need be 
considered. Inspection of (17) and (18) for ¢6=0 shows 


E,(0) = cos? < [Jo(B)+J2(B) |—7 


> ~ 
. p-R; = x sin 8+ ro cos cos 0. 


Then (24) and (25) may be written 


, t/2 Tx S w/2 
E,(6) = sin of cos pike oes ‘ax [ sin y cos y 
—t/2 t =n) 2 
u e 2ikro cos y sin” 812d, (31) 
and 
t/2 Tx w/2 
E(0) = (1 + cos 6) COSH a =je4k> sim Pda: cos? py 
Sati —7/2 
(32) 


-e727kro cos p sin” 912d, 


1953 


Equation (31) involves an integral of an odd function 
from —7/2 to 7/2, and is therefore zero. The first inte- 


gral in (32) is given by 
2t (= sin *) 
cos 
vy 2 


ki sin 6\? 
seni 
Tv 


For the second integral, the substitution D=2kr» sin? 
(8/2) gives an integral which was worked out by Mc- 
Lachlan!® in terms of Bessel and Struve functions! as 


(33) 


TX 
cos. —— elke sin dx Ess 


18 N. W. McLachlan, “Bessel Functions for Engineers,” p. 83, 
Ex. 27, The Clarendon Press, Oxford, England; 1934. 
_ 16 Tabulated in E. Jahnke, and F. Emde, “Tables of Functions,” 
4th ed., Dover Publications, New York: 1945. 
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i} cos? We iP cos ¥qy = “| {40(D) — 
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Then, combining (32), (33), and (34), the normalized 
H-plane pattern given in the text by (9) is 


(1 + cos @) (= sin *) 
COS 
2 v 


ee 2 
oe 
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Ji(D) 
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The Effect of Ions on Magneto-Ionic 
Characteristic Polarization 


W. SNYDERY, 


Summary—tThere is some evidence for the presence of large 
numbers of ions in the E-region. Considering the recent increase in 
the precision with which sky-wave polarization measurements-can 
be made, it seems desirable to fill a gap in the usual magneto-ionic 
theory by investigating the effect on the theoretical characteristic 
polarizations of introducing a mixture of ions and electrons into the 
calculations. A general formula is derived as an extension of the 
magneto-ionic polarization expression of Goubau. 

The 100 ke and 300 kc polarizations predicted by this formula 
have been plotted for a neutral mixture of oxygen ions, nitrogen ions, 
and electrons. The plots show that polarizations are possible in a mix- 
ture that cannot occur when only electrons are present. They also 
show that characteristic polarizations in a mixture do not define the 
medium parameters as completely as they do with electrons only. 
Another interesting feature is that the polarizations in a mixture de- 
pend on the relative numbers of each type of charged particle but 
seem not to depend on the actual charge density when the propor- 
tion of electrons is small. 


rather large concentrations of ions in the E 
layer.4? Usual magneto-ionic theory neglects the 
effect of ions and considers only that of electrons as be- 
ing important. Surely, if there are equal numbers of 
ions and electrons this is a valid procedure, since ions 
are so much heavier than electrons. However, if the ions 
are much more numerous than the electrons they can 
have an appreciable effect on the polarization. 
G. Goubau? has treated the ion-electron mixture ex- 
tensively in terms of the propagation constant. He gives 
an expression for the polarization in a general mixture 
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for zero heavy ion and electron densities, and an expres- 
sion for the polarizations for any ion density when the 
electron density is considerable. He shows that in the 
latter case the electrons chiefly determine the character- 
istic polarizations. 

For the limiting polarization case, he points out that 
the polarizations vary widely depending on the type and 
proportions of ions comprising the mixture. A general 
expression for polarization is not given nor is there any 
discussion of what would happen for any specific choice 
of ions except that if the positive and negative ions were 
of equal charge to mass ratio and equal number, the 
characteristic polarizations would be linear with tilt 
angles of zero and 90° with respect to the component of 
magnetic field in the phase front. 

It seems desirable, now that precise polarization 
measurements are being made, to consider in a little 
more detail what some fairly reasonable mixture of ions 
and electrons will do to the theoretical characteristic 
polarizations. To this end, one should have a general 
expression for polarization in an arbitrary mixture for 
any ionization densities. ; 

Such an expression was derived following the same 
general method indicated by Goubau. That is to say, 
Maxwell’s equations were written for the field with the 


—_ 
conduction current term taken to be >, Nie;V; where 


t 


N, is the number of ith type ions per cubic meter, e; 
charge on ith type ions and V; is velocity of ith type 
ions. A force equation for each type ion was written 


—> 
dV; > 


> >. = 
: =¢(E+ Vi X Bz) — miriVi 


(1) 


Mi 
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where m; is the mass and »; the collision frequency of the 
ith type ion. E is the radio frequency electric field and 
Bz is the earth’s magnetic field flux density. This treat- 
ment neglects restoring forces proportional to particle 


displacement. 
After half a dozen pages of tedious algebra, the ex- 
pression for polarization, Q, was obtained as 


Br?[c? — (1+ a)d| 
2Brc[1 + a+ Bzd| 


Br*[c? — (1+ a)d] \? 
+ = 
/ ear ie ebay oe 


xi(1—72:) ) 
pi Une) 


one 
fj 
th 


Qi,2 ae 


where 


Smt 


pel yeoge— 
— LiNi \ 
Brit yeh oe = 11 — gee) 
Nes 
t= 
Eom wo” 
Bre; 
Sa 
MW 
Vi 
Op oe aa 
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v;=collisional frequency of the 7th type ion 
w=signal angular frequency (27f) 
N;=number of ith type ions per cubic meter 
e;=charge on zth type ion 
m;,=mass of 2th type ion 
Br=flux density of earth’s magnetic field (MKS) 
B,=transverse component of Bz. 
By,=longitudinal component of Bg 
€)= dielectric constant of free space (MKS) 
Q=Hz/Hy where Hy is the magnetic component of 
signal perpendicular to plane determined by By 
and direction of propagation, Hz is the compo- 
nent in plane. In an x, y, z, coordinate system, x 
is the direction of propagation. (See Fig. 1.) 


In the limit x,—0, (1), agrees with Goubau and if only 
electrons are considered it becomes the usual magneto 
ionic polarization expression. 

If we write 


1 “ By*[c? — (1+ a)d] 
2Brc[(1+ a+ Bd] 


Ze gaa ye 
Ova= : 4/ ey 
EAD cay ) pes 8) 


then 
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Equation (3) has been plotted in detail for Q=a+ 76 
with Z and X as parameters and also for ratio of minor 
to major axis versus tilt of major axis of the ellipse rep- 
resented by Q with the same Z and X as parameters.* To 
see the behavior of the ellipse we need only solve for Z 
and X and then, from the plots mentioned, pick off the 
appropriate Q or ellipse. Unlike the case with only one 
type of charged particle, both Z and X are functions of 
both x, and z; and the relationships are such that in the 
general case considered here, no simple normalization 
could be found but, instead, special cases had to be com- 
puted. From these special cases, however, it is felt that 
the effect of ions is fairly clear. 

If we require that the ionosphere be electrically neu- 
tral, we must consider only cases in which there are as 
many negative as positive charges. The most likely situ- 
ation is one having both positive and negative ions, but 
with ions of one sign possessing a higher charge-to-mass 
ratio than the other. A fairly reasonable model might be 
one in which there were free electrons, nitrogen mole- 
cules with one electron missing, and oxygen molecules 
with one electron added. Using Goubau’s mixture co- 
efficient, K, for this mixture we have 

X1 


to +g 


where the subscript 1 refers to electrons, the 2 refers to 
nitrogen molecules and the 3 to oxygen molecules. 


K= 


re 


x 


Fig. 1—Coordinate System (Propagation in +x direction). 


We can obtain some insight into the effect of ions by 
calculating polarizations as a‘function of K keeping 
Xi+X2+x3 constant for particular cases of signal fre- 
quency, magnetic dip angle, collisional frequency, and 
gyro frequency. Z+jX for the cases of f=300 ke, x1 
+x2+2x3=0 and 0.5, 2:=0.27 and 1.06 are tabulated in 
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Ix= 0.5 
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Fig. 2—Polarization of one downcoming magneto-ionic component vs. K at 300 kc looking in the direction 
of propagation in the northern hemisphere. 


TABLE I 
Z4+jX vs. K at 300 kc 


Xi +x. +x3=0. 5 


0.4923-+j1.234 


0.4925 +j1.859 


18.56-+ 784.61 


—54.1—j14,205 


K xi+x.+x3=0 xi+x2+x3=0.5 x1 +22 +23 =0 
z1=1.06 z,=1.06 a=0. 2 =0.27 
1 1.964+ 71.851 1.940-+ 71.208 0.4910+ 71.8506 
1+ 0.491134 71.8514 
5 1.972+ 1.862 1.972+j1.843 0.4929-+ 71.859 
Loe 2.0444 71.959 2.044-+71.957 
108 2.768+ 72.932 2.768+ 72.932 0.6821+ 72.731 
TOs 9.398+ 711.94 9.398-+ 711.94 2.405 + 10.66 
105° 83.58+ 7110.3 83.58+ 7110.3 19.84+ 790.42 
10-7 941.5471, 116.5 922.2+4j1,212 
10-8 —5,833—j21,501 21,220—j9,751 11,088+j18, 794 
0 —191.8—714,013 —3,467 —713 ,953 —49,691—j14,542 


Table I and the corresponding polarization ellipses plot- 
ted in Fig. 2. 

Z+jX for the cases of f= 100 ke, x1+%2+%3=0, 21=0, 
1, and 3.2 are tabulated in Table IJ and the ellipses plot- 
ted in Fig. 3. 

In all cases an electron gyro frequency of 1.3 mcanda 
magnetic dip angle of 62° is assumed. The collisional fre- 
quencies of nitrogen and oxygen ions are taken to be 
v= 1/40r, and v3=1/43n, respectively, following Gou- 
bau. 


Two cases of limiting polarization at 7 mc were calcu- 
lated, one for K =10-* (K drops out of the expression if 
it is greater than 10-%) and the other for K=0. The po- 
larization at 7 mc was circular for the two limits of K so 
no further calculations were made as the trend of the 
curves indicates that in all cases the limiting polariza- 
tions at 7 mc are circular. 

As can be seen from the plots and tables, the trend of 
all the cases plotted is the same. That is, for values of K 
greater than about 10~* the polarization is very close to 
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Fig. 3—Polarization of one downcoming magneto-ionic component vs. K at 100 ke looking in the direction 
of propagation in the northern hemisphere. 
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TABLE II 
Z+jX vs. K at 100 ke 
K X1+x2+x3=0 Xi +X. +x3=0 +o +23 =0 
F 2,:=0 21> Zi ie, 
1 0.616+ 7.616 1.96+j0.616 
10-2 0+ 7.6186 0.619+ 7.619 
10+ 0+7.6433 0.642-+7.646 
10-3 0.875+7.914 2.79+ 71.10 
5(10)-*4 1513-71. 21 
1054 0+ 73.336 3.21+ 73.61 10.23+75.459 
S(10)e° 4.89-+76.46 18.53-++j10.28 
10> 29.4+730.6 86.1+748.1 
10-6 1 293 +7300 1,030 +7372 
TORG 601 —748.2 2,680 —75 ,059 
0 0—j4,810 —62—j4,890 | —193—74,281 


that determined by the electrons only, but as K de- 
creases (proportion of ions increases) the ellipse tends 
to become more circular and the tilt angle changes. For 
very small K the ellipse is essentially circular and al- 
though as K approaches zero the tilt increases, finally 
getting into the “forbidden” quadrant for electrons, the 


ellipse is so nearly circular that the tilt angle is not ob- 
servable. 

One result that comes out of this is that some polari- 
zations, impossible in a medium which has only electrons 
as the charged particles, are possible when ions are pres- 
ent. In fact, all the limiting polarizations (again defined 
here as those corresponding to all x;=0) in a mixture 
could be accounted for only by negative densities of 
electrons alone. Thus, if polarizations were observed 
which would indicate negative electron densities, it 
would not necessarily mean that magneto-ionic theory 
did not apply; it could mean that the medium contained 
a mixture of ions and electrons. aan | 

In addition, one can say from the fact that curves of 
constant collisional frequency for which }°*;=0.5 and 
those for which }>x;=0 essentially overlap for most of 
their length that a given polarization does not uniquely 
determine K or }°>x; in general. a 

On the other hand, this overlap shows that collisional 
frequency could be obtained in principle from a polari- 
zation observation of a magneto-ionic component if con- 
ditions were such that magneto-ionic theory were ap- 
plicable. The procedure for doing it would be cumber- 
some. It would consist of plotting curves for a family of 
collisional frequencies and then interpolating for the 
values of collisional frequency corresponding to the ob- 
served polarization. 
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Table I shows most clearly the extent of the overlap 
for different dx. For the range of K <10~ or so, for 
which the polarizations are not circular, the polariza- 
tions for }°x;=.5 and > \x;=0 seem to be the same if K 
is the same. 

It should be emphasized in concluding, that for the ef- 
fects of ions to be discernible in polarization measure- 
ments, the relative number of ions as compared to elec- 
trons must be enormous. The quantity K may be writen 
for the mixture discussed as 


Ni/m 
> / 

K = — = 
2 N 2/ M2 + N3/ms 


(4) 


where m2/m,= 51,200, m3/m,= 59,200 and since N2=N3 
very closely for values of K considered, we may write 


My ( 1 
muifecets 3 
Neo \mi/m2z + m,/m3 
Mi 
or K = — 3(10)4 
IVe 
Ny Z 
eK = X 6(10)4 since 2V2 = N2+ Ns, (5) 


No + N3 


so that if K=1, N2+N3=6(10)4 Mj. 

It is seen that for this value of K the limiting polari- 
zation is essentially that due to electrons only. In fact, 
the ions do not begin toaffect the polarizations much un- 
til they outnumber the electrons about six million to 
one. This is a much less sensitive quantity for detecting 
ions than is the critical frequency separation. According 
to Goubau, the conditions for reflection in a mixture are 


%1 + X2+ x3 = 1 

X1 
Spite a Seite es 
1+ 


l| 
— 
. 


(6) 


We may write x1 =,"/w”, so that, if w,° is ordinary criti- 
cal frequency in the mixture and w,,” is extraordinary 
critical frequency, 
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Now 
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2” + wy” 
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Hence 
wo 7(1 a 1/K) a Wm”, 
and, in terms of the ordinary critical frequency, 


” 2 
Kam 


1 


Thus 


2 2 Kon” ( Sh 
Gm Ae Way ), 

Bein’ isto ys 
so that, for K =1, the separation is changed by a factor 
of 1/2, for the same ordinary critical frequencies in a 
mixture, in comparison with the pure electron case, even 
though the polarizations are not changed. However, the 
0 and x waves may be reflected from different heights 
where different K’s prevailed. In fact, the reflection is 
likely to occur at a place of few ions, but the polariza- 
tions may be determined in a lower region where ions 
are more likely. The figures given by Seddon? show that 
a ratio of (N2+N3)/Ni=2.8(10)® at least exists at 90 
km, corresponding to K =2(10)7!. K could have been 
much smaller since this value of K is based on J, being 
just below the minimum detectable. 

Since VN, could have been much smaller this shows 
that the effect of ions on the polarizations may have to 
be taken into consideration in interpreting polarization 
experiments. 
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The 1953 Symposium on Tropospheric Wave Propagation 
Within the Horizon at the U.S. Navy 
Electronics Laboratory 
W. C. HOFFMAN* 


Summary—A symposium on ‘‘Tropospheric Wave Propagation 
Within the Horizon” was convened at the U. S. Navy Electronics 
Laboratory during March 30—April 2, 1953, under the chairmanship 
of Dr. Charles R. Burrows. There were three sessions, devoted re- 
spectively to the effects of rough terrain, to the effects of atmospheric 
refraction and to propagation phenomena arising from combination 
of these effects. 


HE SESSION ON terrain effects was about evenly 
Pee between the reporting of experimental re- 

sults and exploring theoretical methods for solving 
the electromagnetic field over a rough boundary. Mr. 
W.S. Ament outlined three possible approaches to the 
latter: (1) the equivalence between wave propagation 
over a rough surface and: propagation through a 
“blobby” medium; (2) approximation of the rough sur- 
face by an infinite set of parallel plates with random 
spacing; and (3) the application of noise statistics. Mr. 
W. C. Hoffman spoke on the approach that has been at- 
tempted at NEL, consisting of solution of an integral 
equation for the surface currents together with applica- 
tion of some recently developed results of probability 
theory. 

Dr. Victor Twersky described his work on reflection 
of electromagnetic waves from a striated surface, using 
the model of an infinite set of randomly spaced semi- 
cylindrical bosses on a perfectly conducting plane. He 
finds that multiple scattering can be neglected with such 
a model at all angles of incidence for horizontal polari- 
zation. For grazing incidence the model behaves like a 
perfect reflector regardless of polarization or the dielec- 
tric properties of the bosses. An interesting result is that 
approximations for reflection coefficients and differen- 
tial scattering cross-sections which conserve energy may 
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be found by retaining only the leading term of an expan- 
sion in terms of 1/k5 cos a (where 0 is the average sepa- — 
ration of neighboring scatterers and a is the angle of in- 
cidence) provided kb cos a>>1. For grazing incidence, 
however, this approximation breaks down and multiple 
scattering must be taken into account. When the radius 
of the cylinders is large with respect to wavelength a 
polar plot of the reflected radiation is roughly semi-cir- 
cular except for a sharp spike-in the specular direction. 
The back scattering cross section in the case of horizon- 
tal polarization behaves essentially the same whether 
the radius is large or small with respect to wavelength. 

Mr. H. W. Smith gave a paper entitled “Power Spec- 
trum Estimates of Radio Waves and Simultaneous Sea 
State for An Overwater Path.” The Tukey method was 
used in estimating the power spectra. There was rea- 
sonable agreement between the power spectra of radio 
field and sea state at centimeter wavelengths but not at 
millimeter wavelengths. A paper by L. G. Trolese and 
J. P. Day, “Diffraction Effects on an Optical Short 
Wave Link,” described the pronounced enhancement by 
night-time refraction of knife-edge diffraction from a 
small hill about a half-mile in front of the receiving an- 
tenna. The possibility of strong local propagation ef- 
fects from portions of an optical path in the neighbor- 
hood of the antenna was also brought out in K. Bulling- 
ton’s paper on “Reflection Coefficients of Irregular Ter- 
rain.” Over one path a change from 16 to 11 db in the 
depth of the null resulted from moving the receiving an- 
tenna laterally a hundred feet. A plot of reflection coef- 
ficient versus roughness parameter for a number of 
paths across the United States showed no particular 
correlation between the two, the mean reflection coeffi- 
cient being 0.28 at 4000 mc. 
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Mr. A. P. Barsis reported on “Tropospheric Propaga- 
tion Measurements within the Radio Horizon—Chey- 
enne Mountain Path.” This is an optical path in Colo- 
rado essentially duplicating air-to-ground propagation. 
A least squares fit of the terrain over that portion of the 
path within line-of-sight of the receiving antenna was 
used to give an effective earth’s radius and smooth- 
earth theory then applied. Rayleigh’s roughness crite- 
rion in terms of the mean deviation of “fitted” from 
“true” terrain height was computed. It was found that 
when the criterion was small, the measured field ex- 
ceeded the computed field; when it was large, the com- 
puted field was the larger. Dr. H. W. Swarm described 
some VHF measurements over extremely rough terrain 
in the Seattle area and a method of computing the field 
on the basis of reflected, refracted and diffracted com- 
ponents. 

Dr. J. B. Smyth led off this session with a comparison 
between a link partially over water from Mt. Wilson to 
San Diego, thus simulating air-to-ground propagation, 
and one wholly over water from San Pedro to San Diego. 
There was no apparent relation between instantaneous 
signal levels on the two links. Mr. J. P. Day then dis- 
cussed a paper titled, “Propagation Characteristics of 
Microwave Optical Links,” concerning fades on links in 
the vicinity of San Diego, San Francisco, and Norfolk. 
On one occasion the height-gain curve for the San Diego 
link under trapping conditions was such that maxima 
under standard conditions became minima and vice 
versa. Thus a receiving antenna located at a normal 
maximum of the height-gain curve would show a deep 
fade. The shift in the height-gain curve could be com- 
puted on the basis of an effective earth’s radius of —2. 
From a continuous record of field strength it was also 
found that there was a fading period of approximately 
14 seconds coinciding closely with the period of ocean 
swells originated by southern hemisphere storms. Of the 
two paths in the San Francisco area, one had a strong 
ground-reflected component and exhibited deeper fades 
than the other. The Norfolk link was notable for the 
number of lengthy “flat bottom” fades observed there, 
some of which were associated with frontal passages. 

Dr. L. H. Doherty gave a paper on “The Effect of 
Atmospheric Ducts on Line-of-Sight Transmission” in 
which he developed ray theory for optical transmission 
for various types of ducts and various transmitter 
heights and proposed supplementing ray-tracing by use 
of the Airy integral in the neighborhood of caustics. 
Mr. Wong described the work at the Aircraft Radiation 
Laboratory in ray-tracing in a paper, “Correlation of 


Radio Holes and Fading Regions with Calculations.” 


Ray tracing with a REAC gives a graphic display of 


-radio-holes for a number of transmitter heights and a 


stratified atmosphere or linear refractive index profile. 


The location of the radio-holes is quite sensitive to trans- 
_ mitter height. It is hoped to employ Freehafer’s method 


and the REAC in a more general wave treatment. 


eye 
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There were a number of comments on the ray-tracing 
approach, which may perhaps best be summed up in the 
statement that the method has practical value but does 
need to be supplemented in regions near cusps and caus- 
tics. To have complete correspondence between theory 
and practice it is necessary to know much more than we 
do about the meteorology of progressive change in the 
refractive index profile in time and space and, in some 
situations, to take into account the ground-reflected 
ray. 

The session on refraction was marked by a lively con- 
troversy between the advocates of ray-tracing methods 
for determination of “radio holes,” on the one hand, and 
the proponents of multipath and focusing-defocusing 
effects as explanations of fading on optical links, on the 
other. The matter of eventually explaining and/or pre- 
dicting fades on optical communication links (and ra- 
dio-holes in the case when one or both terminals are air- 
borne) was probably the main goal of the session. 

Mr. Bullington then described the experiments of 
Crawford and Jakes (BSTJ; January, 1952) on fading 
on optical links caused by multipath transmission and 
showed some highly interesting films of the behavior of 
the received pulses. There was a lively discussion of the 
explanation of the observed phenomena. 

Mr. P. Chavance then described his observations on 
a grazing microwave link in the north of France. The 
median signal fell between the values for diffraction by 
a spherical earth and diffraction by a knife edge. There 
were marked diurnal and seasonal trends in the fre- 
quency of deep fades. It was notable that interference 
type fades could be explained on the basis of two rays 
of nearly equal amplitude and a third ray of much 
smaller amplitude, whose relative phases are slowly 
varying linear functions of time. Mr. Chavance also de- 
scribed a new link consisting of two convergent paths 
over the Mediterranean, from which he expects to have 
some interesting data shortly. Dr. C. C. Aurell spoke on 
fading observed at 3.5 cm on an optical overwater path 
in the Baltic. Space diversity reception was employed 
with one receiver situated at a height-gain maximum for 
a 4/3 earth and the other at a minimum. Two different 
types of fades were observed. The first consisted of deep 
fades on both receivers with no apparent correlation, the 
second, of low signal on one receiver but strong recep- 
tion on the other. The first type may be attributed to 
convergence of the reflected ray (so that its amplitude is 
approximately that of the direct ray) together with a 
third component reflected from an elevated layer. He 
also has noted that deep fades on overwater paths are 
often associated with frontal conditions. 

It seemed to be the concensus that the variation in 
optical fields is due to either multipath transmission or 
focusing-defocusing, and it was suggested by Dr. Smyth 
that the focusing-defocusing effect would appear to have 
a different frequency dependence than the multipath, so 
that this might provide a means of distinguishing ex- 
perimentally between causes of fades. 
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Mr. W. E. Gordon gave a paper on “Radio Scattering 
on Optical Paths,” which is an extension of the Booker- 
Gordon theory of turbulent scattering to the case of 
large blobs (of the order of 10-50 meters) and short 
paths (30-60 kilometers). These new features necessi- 
tate taking account of the non-uniform illumination of 
the blobs and of the searchlight portion of the scattering 
polar diagram. The theory is not intended to apply to 
deep fades but only to the small scintillations on an oth- 
erwise steady signal and, as such, points a way to deter- 
mining the scale of turbulence experimentally through 
measurements of the three parameters: fading ampli- 
tude, fading rate and correlation distance. 

Mr. A. P. Barsis described some measurements on 
fading over two optical links in Eastern Colorado, the 
fades being of the flat-bottom type. His observations 
also exhibited pronounced diurnal and seasonal trends 
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in frequency of occurrence, and the fades were much 
more frequent and severe on 1046 mc than on 100 or 
192.8 mc. Mr. Doherty showed how flat-bottom fades 
might be due to a radio hole. An alternative explanation 
for the fades observed by Mr. Barsis was suggested by 
Mr. L. J. Anderson. Night-time drainage of cold air 
down a valley in front of the farthest receiver could pro- 
vide a reflecting surface and thus give a reflected com- 
ponent annulling the direct ray. For the nearer receiver 
the local terrain was such that a substandard condition 
could be produced by moisture condensation and it was 
pointed out that such a condition would have a much 
greater effect at the higher frequency, in accord with the 
observations. A final paper, “Summary of A.P.O Micro- 
wave Propagation Measurements on Sydney-Melbourne 
Route,” by Mr. J. H. Reen was read by title, describing 
some fades apparently due to substandard conditions 
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The IRE Professional Group on Antennas and Propagation is grateful for 
the assistance given by the firms listed below, and invites application for 
Institutional Listing from other firms interested in the field of Antennas and 


Propagation. 


THE GABRIEL LABORATORIES, 135 Crescent Road, Needham Heights, Massachusetts 
Research and Design of Antenna Equipment for the Workshop Assoc. and Ward Products Div. of the Gabriel Co. 


POLYTECHNIC RESEARCH AND DEVELOPMENT COMPANY, INC., 55 Johnson Street, Brooklyn 1, New ° 


Microwave Precision Test Equipment—Design, Development, Production. 


WHEELER LABORATORIES, INC., 122 Cutter Mill Road, Great Neck, New York 


Consulting Services, Research and Development, Microwave Antennas and Waveguide Components. 


The charge for an Institutional Listing is $25.00 per issue or $75.00 for four 
consecutive issues. Application for listing may be made to the Technical 


Secretary, The Institute of Radio Engineers, 1 East 79th Street, New York 
21, New York. 


